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Abstract. We consider the generation of sound waves in the air medium between 
electrodes at the voltages near electrical breakdown in the presence of the time 
dependent constituent of the electric field. Within the standard multicomponent 
^\f , hydrodynamic model of the weakly ionized gas it is shown that the generation of 

2 I sound is possible due to instantaneous character of the ionization equilibrium. The 

influence of the electronegative ions on the sound intensity is also discussed. 
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1. Introduction 

The acoustic devices based on the plasma and the corona discharge have important 
advantage of direct transformation of the electric signal into sound. This opens up 
possibilities for designing the "membrane-free" acoustic systems [H El [3] in contrast to 
the standard loud-speakers |4j. The main principle of operation of the electroacoustic 
transducers is the transformation of the electric energy into oscillation of the air 
density (pressure). Such transformation can be either direct or indirect. The direct 
transformation is when the electric energy transforms into the motion of macroscopic 
volumes of the media producing the change of the pressure. The indirect one is connected 
with the change of the internal energy, i.e. the heat dissipation by the Joule losses and 
the corresponding thermodynamic process of changing the pressure (density). With 
respect to these mechanisms the transducers are divided into two classes: the "cold" 
and the "hot " ones [5], |6] . The heating mechanism of the sound generation is typical for 
the high frequency corona type and plasma sources [21 [71 E]- The force mechanism of 
the sound generation is characteristic for the low frequency (< 10 kHz) corona sources 
[H El [8]. The simplified electric circuit model was developed to describe the main 
characteristics of the transducers [H |9]. Basically monopolar and dipolar directivities 
are associated to the ionization and drift regions, respectively. The physics behind the 
above-mentioned devices is known and the main predictions concerning the directivity 
pattern match very well with the experimental results. In particular, the effective 
electrical circuit models were proposed [U |7]. They give a clear physical picture on 
the phenomenological level and describe the main characteristics rather good using the 
parameters of the effective electric circuit. Nevertheless, there is the need for the more 
detailed description of the sound generation based on the more detailed consideration. In 
particular, the phenomenological approach fails to answer the questions about the ozone 
generation, the dependence of the characteristics on the geometry and the chemical 
composition. 

The possibility of the monopolar source for the cold corona type transducers 
which would be indistinguishable from the heating related field has been noted in 
[7]. The glow discharge SLT-Ionophone designed by F. Fransson [TOl [11] seems to 
display such possibility. Besides, the fact that the discharge is insensitive to the 
magnetic field signifies that the ion wind effect does not matter for the sound generation. 
Unfortunately, the specific physical mechanism behind this device was not given by the 
authors of [TOl[TT]. Obviously, the "corona wind" field source is incompatible with what 
is known for this dc glow transducer [9]. 

In this paper the theoretical scheme of calculation of main characteristics of the 
sound wave excitation process in the gap between electrodes in the presence of the 
alternating electric field constituent is proposed. We use the standard three-component 
hydrodynamic model of the weakly ionized gas, demonstrating that the generation 
of the bulk charge is the monopole sound source which is generally present in cold- 
plasma transducers. Therefore we identify the physical process of the monopole sound 
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generation in the force driving corona or glow discharge. In particular, this gives the 
explanation of the action of SLT-Ionophone [11]. We also propose our own transducer 
which is based on this effect. 

The paper is structured in the following way. In Section 2 we derive a set of 
equations describing the electroacoustic system under interest, including electrodynamic 
and kinetic equations. Next Sections 3 and 4 are devoted to the experimental setup as 
well as characteristic values and some simple estimates of main quantities appearing 
in the problem, while in Section 5 we discuss the volt-ampere characteristic and 
autoelectronic emission in a concrete device. Our main results are summarized in 
Conclusion. Finally, in Appendix we turn to the simplest possible one-dimensional 
case and obtain a numerical solution of the problem. 

2. Set of equations for electroacoustic system 

The description of the system of interest is based on the system of equations which 
incorporate all main physical processes. They describe the distribution of the electric 
field and charges. This set of equations is well known and widely used in modeling of 
such devices [121 [131 [H] . 

Because of weak ionization of the air medium, we shall calculate different 
quantities, characterizing the electroacoustic system, within the limits of quasi- 
stationary approximation and the corresponding electrodynamic problem. For the 
medium with the conductivity a we have the following system of equations: 

Aip = p, 

| + Vj = 0. (2.1) 

J = je+j-+j+ = OrE, 

where A denotes the Laplace operator; cp and E = —Vip are the potential and the 
strength of the electric field respectively; e denotes the relative permittivity; the current 
density j consists of three terms je, j- and j+, corresponding to electrons, negative ions 
and positive ions respectively; finally, the charge density 

p = e(n_|_ — ?T,_ — rie) , (2-2) 

where e denotes the elementary charge; n+, n_ and rig are densities of positive ions, 
negative ions and electrons respectively. Here and in what follows we use the CGS 
system of units. In particular, e ~ 4.8 x 10~^° un. CGS. 

The above-mentioned quasi- stationary approximation means that in the zeroth 
order the strength E = Eg, where the stationary constituent Eq does not depend on 
time t, while in the first order the strength can be presented in the form 

E = Eo + 6E{t) , (2.3) 

where SE{t) is the alternating constituent (hereinafter also called a "signal") and it is 
assumed that 5E(t) -C Eq. 
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In such a problem setting the main quantity, which is determined by kinetics of 
coUisions, is the effective conductivity a. It is clear that in view of small mobility of 
ions the current is caused mainly by the electronic component. For the same reason 
the ionic component makes the main contribution to the sound pressure. At that the 
substantial factor is the generation of electronegative ions of oxygen O^ in consequence 
of attachment of electrons, that must be taken into account in balance kinetic equations 
and corresponding cross-sections of reactions of birth and annihilation of particles. This 
part of the problem is the standard problem of kinetics and supposes consideration of 
processes of loss and birth of electrons, colliding with molecules and ions [15] . 

To describe the sound generation, one should supplement the electrodynamic Eqs. 
(f l2.ip ) with the equation of motion of the gas (air), taking into account the electric 
force, acting on charges: 

dv VP p 

+ (vV)v = -^— + ^E+..., (2.4) 

dt p p 

where p is the mass density; . . . denotes viscosity terms, which can be neglected under 

consideration of propagation of sound of typical wave lengths 1 -=- 10 cm. Since the 

characteristic frequencies are less than 10 kHz we neglect the heating mechanism of 

the sound generation which is typical for the plasma arc and corona-type devices of 

high-frequency [21 E]. Only the so called force source of the sound is taken into account. 

Linearizing Eq. (12. 4p in a standard text-book manner (see, for example, [IS]), we 

obtain: 

^._Zf,PE^-^.±j. (2.5) 

dt po p po HP 

Here we used the definition of mobility /i: j = /upE, a = up. 

Let us note that p, j, a, n±, Ue and v (as well as Aip, VE and VP) are the quantities 
of the same (first) order of smallness. Similarly the mass density p can be decomposed 
in the following way: 

p = PQ + 6p^ (2.6) 

where po and 6p are quantities of zeroth and first orders of smallness respectively. 
Linearizing the continuity equation 

^ + Po(Vv) = 0, (2.7) 

after some straightforward computations from Eq. (12. 5p we get: 

where c is the speed of sound and jo is the current value in the stationary state. 
We take into account that the mobility depends on the strength of the electric field 
mainly because of the electronic component. Indeed, the electronic mobility /Xg depends 
essentially on E in the interval 1 < P < 5 CGS while the ionic mobilities can be 
considered as constants (see [H]). Obtained Eq. ( (I2.8P ) clearly demonstrates that 
there are three sources of sound of different physical nature. Eq. (12.81) is the standard 
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equation of acoustic oscillations in the presence of the external source (representing 
oscillations of the bulk charge, caused by the presence of the alternating constituent of 
the electric field) . Note that the terms in the right-hand side of Eq. f l2.8p differ by their 
symmetry. The first one represents an isotropic source (/ = harmonics), while two 
others obviously represent the oscillation on the ion wind background and indeed give 
raise to the common dipolar contribution (/ = 1 harmonics) into the acoustic directional 
pattern [Il[7l|8]. The relative magnitudes of these contributions depend on geometry of 
the device. Thus, it is possible that the sound is produced mainly due to the oscillation 
of the bulk charge (the first term in the right-hand side of Eq. (12. 8p ) . Yet it is explicitly 
monopolar in contrast to the statement of ^ that all force sources are dipolar since they 
are connected with the electric wind. From the physical point of view this term is caused 
by the low inertia of the ionization equilibrium. This effect becomes more pronounced 
in the electronegative media where the formation of large amount of anions becomes 
possible. Obviously, the first term exists even if the electric wind is negligible. To our 
knowledge, all "cold" corona sound sources of point-to-plane geometry considered in the 
literature are of dipole nature [HI [1], [71 E] ■ Note that in "hot plasma" loudspeakers the 
monopolar source is due to the thermal heating mechanism [8] . Thus, our analysis shows 
that the monopolar field source term in generally present in "cold plasma" loudspeakers. 
Obviously, Eq. (12. ip and Eq. (14.11) should be supplemented with the equations 
of chemical kinetics, reflecting processes of generation and loss of charges. Under 
the ambient conditions and the electric field strengths of the order of the breakdown 
strengths one can neglect the processes of diffusion transfer and take into account only 
collisional contributions. Thus, the system of equations of coUisional balance in the 
"minimal" model takes the following form: 

' V(r2+/i+E) = l^Y^elT'e - l^+'efl+ne , 



,{«)„ 



dt --^^^^-^ 

Q'fl_ 

-^ - V(n_/._E) = z.^>e , (2.9) 

V(ne/ieE) = \v\\^ - i^'qJJ n^ - v\'^^n+ne 



dt 



"e 1 



where we have taken into account that jj = q-iUi^i and v^ = ±yUjE [i = ±,e); u^ 



(k) 



«j 



are rates of ionization {k = i), recombination {k = r) and attachment {k = a). In 
general, mobilities /ij depend on i? = |E|. The balance Eqs. (12. 9p agree with the 
charge conservation law and in our approximation n^ and Ug are quantities of the first 
order of smallness, therefore the term — z/^gn+ne in right-hand sides is of the second 
order of smallness and, consequently, can be neglected. Eqs. (12. ip . (12. 9p describe the 
sound generation process under interest. This system should be supplemented with the 
boundary conditions on surfaces of electrodes. Note that the resulting sound pressure 
depends nonlinearly on the input signal due to dependence of fi on the strength of the 
field E. The mobility fi itself is determined with respect to the static field, therefore, 
only the static constituent is taken into account in /i(-E). The calculations show (see 
Appendix) that for amplitudes of the alternating field E/ Eq < 0.25 the nonlinearities 
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Figure 1. Scheme of the device. 



are small. 



3. The experimental setup 

The ionic acoustic system (see Fig. [1]) consists of two arrays of thin needles (4) of stainless 
steel of the diameter 0,025 mm with the sequence period compared with their diameter 
attached to two electrodes (2,3) of length 200 mm. Constant voltage (10-30 kV) under 
which the corona discharge exists is applied on the electrodes. While modulating the 
value of this voltage within the limits of steady existence of the corona discharge, the 
generation of the corresponding acoustic signal arises in the interelectrode interval (1). 
Electric measurements of VAC were carried out with the help of the oscillograph LeCroy 
WR-62Xi. The constructed theoretical model describes the really created models of the 
so-called ionic acoustic systems. While working, the improved design, representing two 
electrodes (2,3) in the form of arrays of thin wires of stainless steel of the diameter 
0,025 mm (4), with the sequence period compared with their diameter, was elaborated. 
The series constant-voltage sources (6) and grids (5) were used as sources. On Fig. [2] 
we depict the directional pattern of the device which shows that the first term really 
gives the main contribution at least for the frequencies less than 10 kHz. Measurements 
were carried out in the anechoic box with the help of the gauged electrical-type acoustic 
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Figure 2. Directional pattern. 




Figure 3. Tlie measurement setup of ttie directivity pattern. 



equipment of "Bruel & Kjaer" (see Fig. [3]). The ionic acoustic system was set on 
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the rotary basis, which slowly continuously rotated. Simultaneously the value of the 
amplitude of the generated sound was registered for the frequencies 2, 5, 10 kHz. The 
microphone was located on the acoustic axis of the loud-speaker assemblage, at the 
distance 1 m. The directional pattern has some distortions for the frequency 10 kHz, 
which are explained by commensurability of the wave length and linear dimensions of 
the loud-speaker elements, including the inner grid and the dielectric holders. 

4. Characteristic values and estimates of main quantities 

In this Section we derive some general inferences from the basic equation Eq. ( 12.81) 
and give the estimates for the relevant physical parameters of the setup. It should be 
mentioned beforehand that the characteristic space scales are as follows: L = 2.1cm is 
the distance between electrodes and S = ix 10~^ cm'^ is the typical value of the cathode 
surface area. 

In order to get some estimates we neglect the coordinate dependence of the mobility 
fi and therefore the third term is omitted. This is based on the following reasoning. The 
spatial variation of fi is caused by its dependence on the strength of the electric field 
E. The main contribution to the mobility /i is due to the electronic component with 
the mobility /Xg- For the characteristic values E ~ 10^ V/cm, Pq = 1 aim, under which 
the device is operating, it does not depend essentially on E according to pTT]. This 
justifies the fact that the spatial variation of /i can be neglected. The second term can 
be also neglected provided that jo is small, i.e. no substantial ion wind produced. Our 
experimental setup confirms the possibility of such simplification. 

Bearing in mind these simplifications and using the Poisson equation Eq. ( 12. ip (we 
put £ = 1 for simplicity), from Eq. ( 12. 8 P we get another equivalent form: 

c2 dt^ Anfi dt ^ ' 

This form makes clear another direct inference that the phase shift between the excess 
sound pressure and the voltage of the harmonic signal is equal to — 7r/2. 

The sound intensity can be estimated from the dimension reasoning. Indeed, the 
excess pressure produced by the bulk charge oscillation in the interval of length L 
between electrodes, reads as follows: 

Ps ~ LpEi ~ EqEi , (4.2) 

where Ei is the amplitude of the signal. For the typical values Eq ~ 10 un. CGS, 
El ~ O.li^o we obtain Pg ~ 10 un. CGS = 1 Pa, resulting in the sound pressure level 

This simple estimate shows that in the considered electroacoustic system the sound 
pressure level can be made big enough. Below we refine the estimate basing of the 
complete set of equations. 
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The obtained equations allow to estimate the value of -^ directly from Eq. f l2.9p . 
Using the data of [19], for the characteristic strength of the electric field we get 
z/_^ g ~ 10^ s~^, He ~ 10^ cm^^ and therefore 

^r^bun.CGS. (4.4) 

dt ^ ' 

This leads to the value of the excess pressure Pg — 0.4 Pa and the corresponding sound 

intensity Is = 86 Db. Comparison of the estimates Eq. f l4.2p and Eq. (14. 4p shows the 

self-consistency of the approach. As a result, we have shown that the generation of 

sound is due to the oscillation of the bulk charge density. 

Let us also note that the presence of attachment with formation of anions leads to 

the larger change of the charge density in comparison with the case of the presence 

only of standard ionization/recombination with formation only of cations. Taking 

into account the large sensitivity of cross-sections of anion-aided reactions to the field 

strength, one can draw a conclusion that the presence of electronegative ions promotes 

the power increase of electroacoustic generation (see also [20], where the mechanism of 

intensification of acoustic oscillations at the expense of friction between electrons and 

neutral particles is proposed). 

5. Volt-ampere characteristic 

The primary generation of charge is induced by the static field with the near breakdown 
strength. In order to determine the mechanism of the primary injection of the charges, 
we perform the measurements of the volt-ampere characteristic (VAC) for the stationary 
(no-operation) mode where the sound generation is absent. The VAC of such mode is 
directly measurable. 

The results of the VAC measuring correlate well with the VAC of autoelectronic 
emission (that is explicitly demonstrated on Fig. H]). 

We use the standard expression for the autoelectronic emission current density (see, 

e.g., mm) 

j=j^{Eo/E,Yexp(-^] . (5.1) 

In reality j* and E'* differ very much from the theoretical values because of influence of 
the surface roughening on the work function [22]. In general these parameters strongly 
depend on the geometry of the electrodes. Integrating Eq. (15. ip over the surface of the 
electrode, we get the VAC for the total current: 

UU) = I.[ffe.,[-^), (5.2) 

where h and f/* are some fitting parameters which depend on the distance d between 
the electrodes. It should be noted that the demonstrated correlation between the 
experimental VAC and the theoretical autoelectronic emission VAC Eq. (15. 2p allows 
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Figure 4. Comparison of results of approximation of the experimental VAC by 
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Figure 5. Dependence of parameters [7* and /, on d. 



to use it as a boundary condition in the solution of the problem (which is discussed in 
detail by the example of the one-dimensional model in Appendix). 

With the help of the experimental data on VAC we are able to analyze the 
dependence of parameters /*,f/* on d. The results are shown on Fig. |5l Since both 
fitting parameters J* and f/* depend on distance d there is the relation between them 
which is determined by the specific geometry of the electrodes. With sufficient accuracy 
the following relation holds: 



h{d)/h{d^i^)^0.33 



(5.3) 



with dr, 



9 mm and d^ 



35 mm. It is shown on Fig. [61 The dependence 
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Figure 6. Comparison of [/» and tlie approximation Eq. (|5.3|) . 



h{d) ~ d~^ can be expected from the following reasoning. It is well known that 
for the needle-plate geometry the current depends on the distance as / ~ d~'^ [T7] . 
Obviously, when the roughening is essential for both electrodes, the current depends on 
d as / ~ (i~^, which is in agreement with (f l5.3p ). 

Conclusion 

In this paper we have described the sound wave generation in the gap between electrodes 
in the presence of the alternating electric field constituent. The consideration is based 
on the simple three-component hydrodynamic model of the weakly ionized gas. It has 
been shown that for the cold plasma sound emitters there is the oscillation of the bulk 
charge monopolar source of the sound. We have demonstrated that the bulk charge 
serves as the sound source and constructed the theoretical scheme of calculation of the 
sound and its main characteristic: the excess sound pressure. To test the approach we 
give the solution for the ID model. The theoretical model proposed in the paper is 
aimed at the description of "membrane- free " acoustic systems which have no essential 
disadvantages inherent to standard electrodjTiamic loud-speakers. 



Appendix A. One-dimensional model 

The case of the linear geometry is the simplest one, since all quantities depend only on 
one spatial coordinate (for example, x). Here we consider ID analog of the theoretical 
model described above and calculate the corresponding acoustic pressure. Despite the 
quite artificial restriction of planar ID geometry it could give useful results as shown in 
[H] where the idea of effective thickness and the corresponding reduction of 2D problem 
to ID one has been proposed. Yet we do not resort to the effective thickness approach 
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since the detailed theoretical model is out of the scope of the present paper. 
In the ID case from Eq. ^1^ and Eqs. ((O), ((gl])) we obtain: 

— — = 47re(n+ - n_ - ne) , (A.l) 

1 d'^P d'^P e d 



and 



c^ dt"^ dx"^ /i dt 
dn 



- — (n+-n_ -ne), (A.2) 



dt 
dn_ 

dn^. 



+ Q^ (n+fi+Ex) = 


= ^+,ene - 


(r) 

- i^+;en+ne 






ia) 






[UelieEx] = 


= (41 - 




(r) 



dt dx^'^^"^^"^' V^'' ^u,ej-. ^^.e-^-r--^., (A.3) 

respectively. If the distance between the electrodes is equal to L = 2.1cm and the 
voltage amounts to ipo = 20x 10^/300un. CGS (or 20 kV), then the characteristic value 
of the strength of the electric field amounts to approximately 32 un. CGS (or 9.5 kV/cm). 
According to [H], //+ = /i_ = 1.5 x 10"^ x (3 x 10^) = 450 wn. CGS, 

fie ^ ( 0.027 + ^ °'^^^ ) X (3 X 10^) un. CGS , (A.4) 



V v/l + 0.3^0/ 

where Eq denotes the absolute value of the continuous electric field strength being also 
expressed in CGS units, z/J'i ^ lO^^^ ^ iQe = IQ-^, 



uX'e ~ 8 X 10'^ X exp ) s 

SEq 

»^o....n9_../ 7.32xl0^\ . 



z/^ j ^ 9.54 X 10^ X exp — '■ s'^ . (A.5) 

3-En 



For the electric field strength characteristic value 32 un. CGS we get /Xg ~ 1.85 x 10^, 
i^+]e ~ 8.5 X 10^ and z/Jj ^ 4.6 x 10^. 

Using the results of VAC measurements, for the value 5 = 4 x 10~^ cm"^ of the 
cathode surface area we solve Eqs. ( flA.ip ) and ( flA.3P ) numerically for stationary mode. 
The spatial distributions of species and the total charge are in Figs. IA1IIA2[ 

The non-stationary mode is realized by applying the voltage of the form 

(p = ipQ{l + 0.4 sin(u;t)) , (A.6) 

where u) is the signal frequency, which for illustrative purposes is assumed to equal 27r s^^ 
(so, the period T = 2Tr/u = Is). We solve Eqs. ( (IA.ip )-( flA.3l) ) on the time interval 
[0,2T]. The 3D distributions of charges in {x,t) variables are shown in Fig JA3|A4[ 
According to the results, the corresponding sound intensity is 

/, = 201gf^ (A.7) 
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Figure A2. Densities of electrons (a) and total charge (b) in stationary mode (see 
text). 



where the pressure can be estimated as 
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The calculation gives the result Ig ~ 25 dB. 
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